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We show theoretically that a domain-wall annihilation in two-component Bose-Einstein con- 
densates causes tachyon condensation accompanied by spontaneous symmetry breaking in a two- 
dimensional subspace. Three-dimensional vortex formation from domain- wall annihilations is con- 
sidered a kink formation in subspace. Numerical experiments reveal that the subspatial dynamics 
obey the dynamic scaling law of phase ordering kinetics. This model is experimentally feasible 
and provides insights into how the extra dimensions influence subspatial phase transition in higher- 
dimensional space. 



A tachyon is a hypothetical supcrluminal particle that 
violates causality in special relativity. Most physicists 
deny its existence because it is inconsistent with the 
known laws of physics. However, in quantum field the- 
ories, a tachyon field can exist due to the instability of 
quantum vacuum. Here, the 'instability' means that the 
state is at a maximum of an effective potential V(T) for 
the tachyon field T [Fig. HJa)]. The tachyon field grows 
exponentially with time and rolls down toward a mini- 
mum of the potential as the true vacuum. This process 
is called tachyon condensation [l[. 

Tachyon condensation is a key concept used for de- 
scribing the dynamics of string theory, a promising can- 
didate for a 'theory of everything' that describes all fun- 
damental forces and forms of matter in nature 0- In 
this theory, a tachyon exists in a system containing a 
Dirichlet(D)-brane and an anti-D-brane, where the for- 
mer is an extended solitonic object Q and the latter its 
anti-object. The two annihilate in a collision similarly to 
particles-antiparticles annihilation in a collision. The an- 
nihilation process is described by tachyon condensation, 
and the system falls into the true vacuum after complete 
annihilation [l[. 

A remarkable application of tachyon condensation is 
in brane cosmology m which the Big Bang is hy- 

pothesized to occur as a result of a collision of a brane 
and an anti-brane. After this collision, lower-dimensional 
branes remain as relics of tachyon condensation [l[ , which 
correspond to cosmic strings in brane cosmology 0^] 
[see Fig. QXb)]. This situation resembles conventional 
phase transitions accompanied by spontaneous symmetry 
breaking (SSB), resulting in the formation of topological 
defects via the Kibble- Zurek mechanism EH- This 
mechanism produces topological defects in the early Uni- 
verse due to phase transitions [l2| , which has been tested 



in several condensed matter systems I3"21|. In contrast, 
tachyon condensation as a SSB phenomenon has not yet 
been well understood. Because it may lead to defect nu- 
clcation in a restricted lower-dimensional subspace, the 
dynamics should be affected by the degree of freedom 
associated with the extra dimension. However, the in- 
fluence of the extra dimension has never been discussed, 
partly because such phenomena are absent in actual sys- 
tems. 

Here, we provide a groundbreaking system to tackle 
this problem, using atomic Bose-Einstein condensates 
(BECs) . Tachyon condensation is simulated by consider- 
ing vortex formations from a pair-annihilation of domain 
walls, i.e. branes, in binary BECs. This system is ad- 
vantageous in that we can theoretically and experimen- 
tally address the nonlinear dynamics of branes, such as 
defect nuclcation and subsequent dynamics, which is dif- 
ficult in string theory. Anderson et al. [22I ] observed the 
creation of vortex rings via the dynamic (snake) insta- 
bility of a dark soliton in two-component BECs, where 
the nodal plane of a dark soliton in one component was 
filled with the other component, and then, the filling 
component was selectively removed using a resonant laser 
beam. Recently, we proposed that domain walls in phase- 
separated two-component BECs correspond to D-brancs 
in the sense that vortex lines (strings) can terminate 



on them [23|. Therefore, the experiment in Ref. [22 1 
may be interpreted as the demonstration of defect for- 
mations via brane annihilation [S2J , although those phe- 



nomena have never been understood as SSB phenomena 
in a restricted lower-dimensional subspace. In this Let- 
ter, we theoretically show that a domain wall annihilation 
causes spontaneous Z2 symmetry breaking and phase or- 
dering dynamics in the two-dimensional subspace, where 
a tachyon field is introduced by projecting the original 
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FIG. 1: (a) Schematic diagrams of tachyon condensation, (b) Schematic diagrams of an annihilation of a brane pair and the 
relic lower-dimensional objects, (c) Dynamics of a domain wall pair annihilation in a numerical experiment for v = 0.84 and 
A0 = 7T. Panels show the time development of the isosurface = l^l- Black curves represent the cores of vortices in the 
'I'l-component (not shown for t/r = 0). The phase arg^i strongly fluctuates in the region ~ between the branes, where 
there are many vortices with little vorticity of mass current (t/r = 32). The initial fluctuation grows into a meshed structure 
(t/r = 40). The structure yields serpentine vortices, which trap the ^-component along the vortex cores. The ^-component 
propagates along the cores and causes varicose oscillations of the isosurface. Reconnections between vortices accidentally create 
vortex rings or disk-shaped density pulses, which escape to the bulk, (d) Panels show the time development of the projected 
field 4>(x,y,t) of Eq. J5j, calculated from the results of (c). The box sizes are 80£ x 80£ x 25. 6^ for (c) and 102. 4£ x 102. 4£ for 
(d). 



order parameters onto the branes in three-dimensional 
space. Our theory is justified by demonstrating the scal- 
ing law of phase ordering kinetics jS3| in numerical ex- 
periments [see Fig. [Tfd) and Fig. [U[b)]. Although the 
analogue of the brane annihilation was simulated exper- 
imentally on the AB phase boundary of superfluid 3 He 
[26|, its theoretical explanation remains lacking. This 
work provides the first theory of brane-annihilation phe- 
nomena in condensed matter systems. 

We consider two-component BECs, which consist 
of condensations of two distinguishable Bose particles. 
Two-component BECs are well described by two com- 
plex order parameters, &j (j — 1,2), in the Gross- 
Pitaevskii model at zero temperature [Sl[. The order 
parameters in an uniform system obey the action S = 

J dt J d 3 x (ihY^j ^jdt^j - K — Vj with kinetic energy 
density K. = V ■ ^-\V^j\ 2 and potential energy density 



V 



Here, we have ex- 

2Trh 2 ajk/m jk (j, k = 
nj +mZ and the s- 



prcssed the coupling constant gjk 
1,2) with the reduced mass rnj^ = 
wave scattering length ajk between atoms in the tyj- and 
'i'fe-components. The chemical potential fij is introduced 
as the Lagrange multiplier for the conservation of the 
norm Nj = J d 3 x\^j\ 2 = const., which gives the particle 
number of the ^ 3 .-component. We consider strongly seg- 
regated BECs and set the parameters as mi = m-2 = m, 
9u = .922 = g, = v , and 312/3 = 7 = 2 ; this pa- 

rameter setting is experimentally feasible, e.g. Ref. [281 ] - 
The time and length scales of our system are character- 
ized by r = hj Hi and £ = h/ y /mfj, 1 , respectively. 

Let us consider the annihilation of a domain wall 
(brane) at z = —R/2 and an anti-domain wall (anti- 



brane) at z = R/2 perpendicular to the z-axis, be- 
tween which the ^-component is sandwiched by the 
two domains occupied with the ^/i-component. We de- 
fine the inter-brane distance R as the distance between 
the two planes defined by l^il = | \I/ 2 1 - The distance 
R increases with N2, which is controllable experimen- 
tally [22j. Because the 'penetration' of the amplitude 
|^i I (l^l) decays exponentially with distance into the 
vj/ 2 _(vl/ 1 -)domain, the short-range interaction between the 
branes works effectively only when R is comparable to the 
'penetration depth', or the brane thickness, and then, the 
annihilation process can start substantially. The trivial 
process of pair annihilation is that the branes collide to 
leave the trivial state ^1 = const.. However, the an- 
nihilation processes become nontrivial depending on the 
phase difference A6 = arg>I>i(z > R/2) - ai-g^z < 
-R/2). 

Since the essential mechanism of the nontrivial annihi- 
lation has been discussed partly in Ref. (S2[, we explain 
it briefly here [29[. To capture the essence of the anni- 
hilation process, we assume that two branes with large 
i?(^> £) are brought rapidly to a small distance R < £ 
but the two \l/i-domains are disconnected at t = [30|. 
In the annihilation process, junctions connecting the two 
^i-domains emerge in various places on the x-y plane 
around z = 0. For AO ^ 0, a junction causes a su- 
perfluid current along the z-axis with a current velocity 
v ± ~ ^AO/i? > or v± ~ ^(A9 - 2tt)/R < 0. If 
the current velocities through the two neighbouring junc- 
tions are parallel, the annihilation is completed between 
the junctions. On the other hand, the two junctions with 
opposite velocities leave a single-quantum vortex in the 
'I'l-component, where the ^-component is trapped in 
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the vortex cores so that N\ and N 2 are conserved. 

Although the growth rates of junctions with v± > 
and v±_ < are generally different, they are statistically 
equivalent for A© = tt. Then, the junctions grow in 
a random fashion from initial random fluctuations, and 
vortices emerge as serpentine curves along the bound- 
ary between the two opposite junctions. These scenarios 
are demonstrated numerically as shown in Fig. QJc) [H3j. 
The snake instability observed by Anderson et at corre- 
sponds to the coincident limit of the two branes (R — > 0) 
with A© = 7r, where the vortex ring nucleation results 
from the spherical geometry of the external potential [22| ■ 
More generally, we can regard the inter-brane distance R 
and the phase difference AO as two parameters to char- 
acterize the dynamics of the branc annihilation. 

To describe the annihilation process systematically, 
we construct an effective field theory parametrized with 
R and A8. Tachyon condensation in string theory is 
explained by introducing a growing field, that is, the 
tachyon field, in the lower-dimensional space spanned by 
the coordinates along the branes In a similar manner, 
we consider an effective tachyon field T(x,y,t) as a real 
scalar field in a two-dimensional space parametrized by x 
and y. On the basis of the vortex formation mechanism 
explained above, we introduce the variational ansatz, 



tanh 2 (z/a)e ±tAe + Tscch 2 (z/6 



* 2 = yJn 2 (T)sech(z/U), 



,(1) 
(2) 



where T(t = 0) =0 and the sign ± in Eq. (p} takes 
+ (-) sign for z < (z > 0) (HI- This ansatz de- 
scribes the annihilation starting from a small inter-brane 
distance, where the population of the ^-component is 
small between the branes at t = 0. The parameters 
rij = n 2 (t = 0)(> 0) and £j_ are determined to minimize 
the energy J d 3 x(K. + V) at t = 0. The ansatz at t = 
reasonably reduces to the solution of a dark soliton with 
£_l £ and R -> for n§ ->• and A6 = 7r. The growth 
of the tachyon field from T = into T > (T < 0) 
induces a superflow with v± > (v± < 0). Although 
the dynamics may be described more precisely with ad- 
ditional variational parameters, our simplest ansatz is 
enough to capture the essence of the vortex formation 
dynamics. 

The effective potential for the tachyon field T, namely 
the tachyon potential V, is defined as a function of T, 
V(T) = dzQC + V) with d x ^ 3 = dy^j = 0. We as- 
sume that the density parameter n 2 (T) > is determined 
so as to minimize V. It is straightforward to obtain the 
form 



V 



m?6 
9 



(3) 



where F n = A n (A@,i/,i) + B n (AO,v : j)8(n 2 ) with a 
step function, 6(n 2 > 0) = 1 and B(n 2 < 0) = 0. The 



contribution from the gradient V|| = (d x ,d y ) along the 
coordinates parallel to the brane is calculated similarly. 
The original energy is then reduced to the form 



E m = J dxdy [G(T)(eV||T) 2 + V(T)\ 



(4) 



Here, the coefficient G(T) > in the gradient term de- 
pends on T (see Ref. [29(] for details). 

Equation (|4|) represents the effective energy for the 
field T in the 'projected-2D' system. The potential V is 
symmetric V(—T) = V(T) for A6 = tt, because the co- 
efficients Fi and F3 are proportional to cos(A6/2). The 
Z 2 symmetry comes from the degeneracy of the kinetic 
energy oc ~ (hir/mR) 2 for ±T owing to a symmetry in 
the initial configuration; Rc'I' 1 (t = 0) = — Re^ifi = 0) 
with AO = tt. In terms of relativistic quantum field 
theory, we find a particle-like state with a mass mx = 
yJV"{T = 0)/2 = m yJi^F 2 /g by expanding the poten- 
tial V around T = 0. For F 2 > 0, this describes a con- 
ventional particle; however, for F 2 < 0, this describes 
a particle with purely imaginary mass < 0, i.e. a 
tachyon. The existence of a tachyon implies instability 
of the system, because a tachyon rolls down from the po- 
tential maximum at T = toward a potential minimum 
with T > or T < [32j . Figure [Ha) shows the tachyon 
potential V(T) for AO = tt and 7 = 2. The coefficient F 2 
increases with v = /J, 2 /fJ,i, and the tachyon potential is 
convex upward around T = for v < 1. Because R is an 
increasing function of the instability becomes stronger 
as the inter-brane distance R decreases. 

The tachyon field T is an analogue of the order- 
parameter field, e.g. the magnetization density, of a 
ferromagnetic system [33} in a continuum description. 
The rolling tachyon corresponds to the spontaneous Z 2 - 
symmetry breaking from the zero 'magnetization' T = 
toward a macroscopic 'magnetization' T > or T < 
in the ordered phase for AO = 



341. In this sense, 



the inter-brane distance R and the phase difference AO 
(7^ tt) play the roles of the 'temperature' and 'external 
magnetic field', respectively. The coefficient F 2 increases 
with the 'temperature' R, implying that the instability 
becomes weak for a small inter-brane interaction for large 
R. Because the inter-brane interaction decays exponen- 
tially with R for large distance and the instability van- 
ishes precisely for R — > 00, the infinity distance may 
correspond to the 'transition temperature'. Although 
F 2 reached zero for finite R in our effective model, the 
tachyon potential V(T) well reflects the nature of the in- 
stability for small R. On the other hand, the field T feels 
the 'magnetic field' for < AO < tt (tt < AO < 2vr), 
and then, 'magnetization' T > (T < 0) is energetically 
favourable. In the original 3D system, this is energeti- 
cally explained from the difference in the kinetic energy 



induced by the rolling tachyon, 
2 



-Mf) 2 forT>0 



and 



ft(Ae-27r) 



mR 



for T < 0. 
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FIG. 2: (a) Plots of the tachyon potential with AO = tt 
for v = 0.84, 0.92, and 1.00. The potential V(T) takes its 
minimum value V(Tt) at T = ±2],. (b) Structure factor of (j> 
in the numerical experiment [Fig. [TJd)]. The factor S_{q,i) is 
the average of S(q, i) over directions of q. The dashed line 
is the expected law, S(q) cx q~ z for q/ho 2> 1, which is not 
observed in the small system. 



A vortex between junctions with v± > and v± < 
is considered a kink between the regions with T > 
and T < in the projected-2D space. The validity of 
the effective theory is confirmed by evaluating the defect 
nucleation rate. The effective theory neglects the trans- 
fer of particles in the x- and y-dircctions in the defect 
formation process, violating the law of particle number 
conservation. However, the violation influences little on 
the rate evaluation for small R. According to the evalu- 
ations with several analytical and numerical methods in 
Ref. the rate consistently decreases with R. 

Finally, we provide the most important evidence show- 
ing that the field T behaves actually as a two-dimensional 
order parameter. According to the scaling law in the 
phase ordering kinetics 



S3] 



the spatial structure of the 
order parameter is characterized by a single length scale, 
i.e. the mean inter-defect distance, after a rapid quench 
from the disordered phase into the ordered phase. In 
the simulations of the 3D system, the density of kinks 
in the field T is calculated from the line density ^2D of 
the projection of vortices onto the x-y plane by assum- 
ing that there is no overlap between the projection lines 
of different vortices. To visualize the dynamics of T, we 
introduce the projected field 



{x,y,t) 



m 
li 



+00 



dzv± (r, t) 



(5) 



with the superfluid current velocity vj_ = 

l*j| 2 ^ ar S*3 n r i - 1 • ,1 ; 1 

m j2 \^ k \ 2 ■ r i rom kinks m the projected- 

2D space, we have "I^ ^ and then both fields 4> 
and T are constant, <j)(x,y) ~ ±7r and T{x,y) = ±Tj,, 
where V(T) takes its minimum value at T = ±T{,. The 
nonzero "I^-component around the kink cores could 
yield different spatial dependences of the two fields. 
Thus, a spatial structure of the projected field 4>(x,y,t) 
represents that of T(x, y, t) on a length scale larger than 
the kink width (, which is sufficient for the following 



analysis. 

If the field <f> obey the scaling law, its structure factor 
S(q,t) for the wavenumbcr \q\ <C 1/C is written with a 
time-independent function T as ^F{q/l2T>) = S(q,t)l% B . 
Here, the scaling form is expected to follow the universal 
law S(q,t) ~ (^2D/<7) d+1 /'2D' known as the Porod law, 
with the spatial dimension d = 2 for q/l^B ^ 1 [s^ ]. 

Figure HJb) shows scaling plots of S(q,t) for the nu- 
merical experiment of Fig. [TJc). The dynamics of <f> 
in Fig. QJd) resembles ferromagnetic relaxations after 
rapid quenching. In fact, the scaling plots almost co- 
incide with each other after the domain structures of <f> 
become clear. The scaling behaviour is seen from the sim- 
ilarity between the patterns oit/r = 104 and i/r = 176 
in Fig. [TJd). These facts show the domain- wall annihi- 
lation is regarded as phase ordering in the projcctcd-2D 
space. 

The time dependence of Z2D contains a statistical infor- 
mation. The density Z2D follows a power law, 1 /I2D oc t x l 2 
[i^], which indicates that the projected-2D system is dis- 
sipative. The dissipation may come from some degrees 
of freedom, neglected in the effective theory, such as the 
motions of the ^-component along the vortex core, emis- 
sions of vortex rings, and density pulses to the extra di- 
mension induced by vortex reconnection [see Fig. [T] (c) 
and the caption]. It is interesting that strings stretched 
between the brane and the anti-brane cause a complex- 
scalar tachyon field in string theory [l[ while the real- 
scalar field comes from fluctuating fields in the wall-anti- 
wall background in our system. If vortices are stretched 
between the domain walls, a 'vorton', a point-like de- 
fect in 3D, can be nucleated due to superflow of the ^2- 
component along the vortex core 35|. These issues will 
be discussed elsewhere. 

Our proposal gives the first realistic example of non- 
rclativistic tachyon condensation due to the brane anni- 
hilation phenomena. Because the instability is essentially 
caused by the spatially inhomogeneous junctions between 
two domains described with the same order parameter, 
the phenomenon may be generalized to a problem of a 
sudden connection between two media in the same or- 
dered phase. Therefore, it is expected that our theory 
can be extended to various condensed matter systems. 
The 3 He-brane experiment j2(| should be retested from 
the viewpoint proposed in this work. Because techniques 
for detecting the vortex line density in the 3 Hc-system are 
well developed, some statistical information about the 
brane annihilation could be observed from the power-law 
decay of the defect density. 
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SUPPLEMENTAL MATERIAL 
Vortex formations from brane annihilation 

Figure[3]summarizes the brane-anti-brane annihilation 
and resulting vortex formation in our system. The trivial 
process of pair annihilation is that the branes collide to 
leave the trivial state ^1 = const, by pushing away the 
^-component, as shown in Fig. GJa). However, the an- 
nihilation processes become nontrivial depending on the 
phase difference AO = arg^'i(z ^> R/2) — arg^E'i(z <C 
-R/2). 

To capture the essence of the nontrivial domain wall 
annihilation, we consider that the brane and the anti- 
brane are brought close to each other at a small distance 
R at t = 0. In the annihilation process, junctions con- 
necting the two ^1 -domains emerge in various places be- 
tween the branes by pushing out the ^-component in 
the region -R/2 < z < R/2. For A6 ^ 0, the junction 
causes a superfluid current along the z-axis with a current 
velocity v± ~ ^AQ/R > or v± ~ ^(A9 - 2n)/R < 
[see Fig. EJb)]. If the current velocities through the two 
neighbouring junctions are parallel, the annihilation is 
completed between the junctions. On the other hand, 
the two junctions with opposite velocities leave a single- 
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FIG. 3: Schematic diagrams of vortex formation in brane annihilation process, (a), One-dimensional diagram of the order 
parameter amplitudes and l^l- The two ^-domains are separated by the ^-domain (left), and then, a junction is 
formed between the two domains (centre), leading to annihilation (right), (b), Two-dimensional diagram of the two colliding 
branes (planes defined by |$i| = ^2 1). The lower figures represent the phase arg^i(z) and its change of direction between 
2 <C —R/2 and z ^> R/2 in the configuration space. The phase difference AO between the two domains (left) induces a 
superfluid current through the junction with direction v± > (— >) or v± < (<—) (centre). Between the two neighbouring 
junctions, the annihilation finishes or a single-quantum vortex is left (right), (c), Three-dimensional diagram viewed from the 
positive direction of the z-axis. For AO = tt, an initial fluctuation develops into a random distribution of junctions with v± > 
(0) or v± < (®) (left), causing a meshed structure (centre). A serpentine vortex line is formed along a boundary (dotted 
curve) between two junctions with opposite velocity (right). 



quantum vortex in the ^i-component, whose core is filled 
with the ^-component. 

The vortex formation strongly depends on how the 
junction grows from initial fluctuations. Although the 
growth rates of junctions with v± > and v± < are 
generally different, they are statistically equivalent for 
AO = tt. Then, the junctions grow in a random fash- 
ion from initial random fluctuations, developing into the 
meshed structures shown in Fig. EJc). Vortices emerge 
as serpentine curves along the boundary between the two 
opposite junctions. 



Method of numerical simulation 

Here, we briefly explain the methods of our numerical 
experiments. The sequence of the numerical computa- 
tions of time development is as follows: (i) A state with 
a brane-anti-brane pair is obtained by numerically min- 
imizing the energy functional / d 3 x(K. + V) under the 
boundary condition ^1 = at z = with arg'I'i = 
±A0/2 for z ^ 0. (ii) The initial state of the time de- 
velopment is prepared by adding a random noise to the 
stationary solution, (iii) The time development from the 
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initial state is calculated by numerically integrating the 
equations of motion for the action S 1 , called the coupled 
Gross-Pitaevskii equations, 

ft 2 

k 

The numerical integrations are performed under the peri- 
odic boundary condition at x, y = ±51.2£ and the Neu- 
mann boundary condition at z = ±25. 6£. The linear 
stability of a stationary solution <&j(z) of a brane-anti- 
brane pair is investigated by numerically diagonalizing 
the Bogoliubov-de Gennes equations, which are obtained 
by linearizing the Gross-Pitaevskii equations ([6]) with 
respect to a collective excitation S^j(q) = ^i{r,t) — 
= e r " t [u{zy q - r - i "« t - v*{z)e- i '*" r+iM « t ] [Sl|. 

Coefficients of the tachyon potential 



is fixed as £_l = £, for v < §7 + \ since the ^-component 
vanishes with n 2 = 0. 

The tachyon mass mj is written as 

m\ = ~{\ + 2 7 2 + 3u - 6 7 j/)^i. 
1 45 g 

The instability condition m 2 ^ < is satisfied for 1 + 
27 2 + 3v — 67^ > 0. The instability tends to be weak 
as v increases, which is consistent with the numerical 
results of the linear stability analysis. This represents 
that with increasing v the inter-brane distance R ~ £_l 
increases, so that the inter-brane interaction becomes 
weak. Remember that our ansatz is not intended for 
the case of R much larger than the brane thickness. In 
fact, the inter-brane distance R approaches infinity for 
v — > 1 in the numerical results although one obtains 
R = 2arcsinh ^/ gn 2 / ' fiij £1 = \/6arcsinh ^1/3^ £ for 
v = 1 in our analytic model. 



Here we show the coefficients F n of the tachyon po- 



tcntial V = 


J2n-o FnT n for the ansatz used in the 


text. 


The coefficients F n (n = 0, 1, 2, 3, 4) are written as 


F n = 


An + B n 0(n 2 ) with 




A i = 


16 

35' 






2 A6 
- cos , 

3 2 ' 




A 2 = 


— f2 7 + 3cos 2 — -6^-2) , 




A l = 


-i(2 7 -6^ + 7) cos^, 




A = 


^(27-6^ + 7), 




B 4 = 


32 7 2 
75 ' 




B 3 = 


4 2 A9 
— 7 cos , 

5 ' 2 ' 




B 2 = 


3 2 2 A @ 3 2 7 2 32 7 i/ I67 

— 7 2 cos 2 — H ■ -, 

8 ' 2 45 15 45 






Bx = 


-37(27- 6^+ 1) cos— , 




Bo = 


-A(2 7 -6^ + l) 2 , 



where we used 



;i 4 

U^ + 3 



-1/2 



n 2 



Re 



4 7 3 7 ^ A8 2 7 1 

1 1 1 cos h 2u 

5 4 2 3 3 



Surface tension of a kink 

We consider a flat kink parallel to the y-axis in equi- 
librium in the case of A6 = ir with V(-T) = V(T). The 
kink profile is the solution of the equation 

Integrating this equation once, and imposing the bound- 
ary conditions T(x = 0) = and T(±oo) = ±T b , one 
obtains 



IV(T)-V(T b ) 



G{T) 



This result can be used to evaluate the energy per unit 
length of kink, that is the surface tension, as 



— 00 
+00 



dx [G(T)(£d x T) 2 + V(T) - V(T b )] 
dx [V(T) - V(T b )} 



4£ / dT^G{T) [V(T)-V(T b )}. 
10 



The coefficient G is written as 
2 9 //17 2 /32 



G = 



^T + cos^) 8(n 2 ] 
3 64 gn 2 \15 2 / 



fits 
9 



Mi 
9 



The conditions £j_ > 6 an d n 2 = n 2(T = 0) > are 



reduced to ^7 + i < v < ^7 + |. The length parameter 



where we used a function 9(n 2 ) = 1 for n 2 > and 
0{n 2 ) = for n 2 = 0. Then the field T(x) is non-smooth 
at x = ±xo, where n 2 = for |x| > xq > 0. The field 
T should become smooth if the gradient energy of "J 2- 
component is treated appropriately. However, a differ- 
ence due to this unphysical non-smoothness would not 
affect on our order estimation. 
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Defect nucleation rate 




V 1 

0.84 0.86 0.88 0.90 0.92 
v 

FIG. 4: Defect nucleation rates. Vortex nucleation rates 
estimated through different analyses: the analytical estima- 
tion with Eq.((6| (solid curve) and the numerical (squares) and 
analytical (dashed curve) estimations with Eq. (|8]) based on 
the linear stability analysis. These estimations are compared 
with the projected vortex core density I20 (circles) at the same 
time as obtained from the numerical results shown in Fig. [5] 
The inset shows the numerical result of the (/-dependence of 
the inter-brane distance R. 

We can evaluate the vortex nucleation via the brane 
annihilation based on the effective tachyon field theory. 
In the following, we confine ourselves to A6 = 7r. A vor- 
tex 'projected' on the effective 2D system corresponds to 
a kink in the field T. When vortices are localized on the 
2 = plane immediately after vortex nucleation, the to- 
tal length L v of vortices in the original 3D system equals 
that of the kinks in the 2D system. The potential energy 
decreases when the tachyon field grows from T = to- 
ward the potential minima T = ±7],. The total energy 
can be conserved in the process of tachyon condensation 
if the decrease in the potential energy is compensated 
by the increase in the gradient energy by nucleating the 
kinks. Because the energy density in the initial state is 
V(0) in the 2D system, the line density Zkink of the kink 
can be written as 



'-kink — 



7(0) - V(T b ) 



with the 2D system's area A and the tension a (energy 
per unit length) of the kink. Note that the conserva- 
tion of N\ and N 2 is not considered in the estimation 
explicitly. Exactly speaking, in the process of vortex for- 
mation described by the variational ansatz, N\ and N 2 
are typically increased and decreased with time, respec- 
tively A simplest form of the ansatz is introduced so as 



to describe qualitatively the annihilation starting from a 
small inter-branc distance in strongly segregated conden- 
sates. In that case, the changes of N\ and N 2 are small 
in the process described with the ansatz, which causes a 
minor correction to the rate estimation. 

A kink solution of the field T(x) parallel to the y-axis 
satisfies the equation SE 2 d/6T = with the boundary 



condition T - 
calculated as 



±Th for x — > ±00, and the tension is 



dT\J G(T) [V(T)-V(T b )]. 



(7) 



The line density /kink is a decreasing function of v be- 
cause the tension a increases with v more rapidly than 
V(Q) — V(Tb). The increase in the kink tension a in the 
2D system may correspond to the increase in the vor- 
tex tension in the original 3D system, where the occu- 
pancy of the ^-component in the vortex cores of the 
^1 -component is increased. 

To demonstrate the validity of the effective model, 
the analytical estimation above is compared to the nu- 
merical results based on the linear stability analysis of 
the original 3D theory. We investigated the stability 
around the initial state against a fluctuation S^j = 



,r„t 



[u{z 



v*{z)e- iq - r+l ^ t ] in the form of the 



Bogoliubov excitation with growth rate r g , wavenumber 
q, and frequency uj q . The q-dependence of T q was in- 
vestigated in a previous work [S2|. For small q, T q in- 
creases linearly with q and reaches its maximum r max at 
q = g max . A further increase in q leads to decreasing T q , 
which becomes zero at q = qt op - The frequency uj q is zero 
for T q 7^ 0. The values r max , g max , and q top decrease as 
v increases. If the instability is sufficiently strong, i.e. 
r m ax is large, the line density Zkink can be estimated as 

7 Qmax 



kink 



(8) 



The numerical estimation with Eq. ([8]) (squares) is rea- 
sonably consistent with the analytical estimation with 
Eq. ((6|) (solid curve) in Fig. |4] Both estimations are in- 
valid for v — > 1. This is because the numerical results 
show that R diverges and r max approaches zero for v — > 1 
[see the inset of Fig. [4], where both the effective theory 
and the linear stability analysis are inapplicable. 

It is instructive to construct an equation of motion 
for the rolling tachyon from the effective field the- 
ory. The equation of motion conventionally contains a 



term proportional to SE 2 d/ST cx (— G(0)£ V 



(6) around T = 0. By considering a small perturbation 

q for q — > 0, a dimensional 



oc e Tqt sin(gx) with T q 
analysis gives 



Then, one obtains 



-G(0)£ 2 Vjj + m 2 T )T (T«l). (9) 



hT q =tqJj-[G{0)eq 



m 



(10) 
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and 



2G(0)' 



(11) 



which gives a quantitative agreement with the defect nu- 
clcation rate, as shown in Fig. 0|b) (dashed curve). For 
q 0, the equation © yields a Klein-Gordon-type equa- 



VnT, which describes a masslcss particle 



tion d 2 T - 

with an imaginary phase velocity ct = y/grn^/m^i^. 
This particle has a growth rate for zero momenta, T q = 0, 
as a consequence of the energy conservation in the growth 
process. We may construct an equation of motion sim- 
ilarly around the potential minima T = ±Tf,, and then, 
the equation describes a massless particle with a real 
phase velocity. However, the minima T = ±T{, should 
be distinct from the true vacuum of this system, because 
the phase gradient 9 z arg>I'i is localized around z = 0. 
In the true vacuum, *f>i = const, and ^ = 0, obtained 
after the complete annihilation, the dispersion of per- 
turbations has in turn a real phase velocity W [i\jm with 
= V Mi/ 9 f° r the massless particle, that is, a phonon. 



Decay of projected vortex line density 




60 80100 

t/T 



The time dependence of the projected defect density 
?2D also contains information that is important for un- 
derstanding the system. The phase ordering kinetics 
suggests the growth law 1/Izd °c t 1 / 3 or l/hj) °c t 1 / 2 
depending on whether the order parameter field is con- 
served or not, respectively [S3|. The decay of the density 
?2D follows a power law, l/hu ck t 1 / 2 , in a manner similar 
to the different values of v as shown in Fig. [SJ This in- 
dicates that the effective energy in the projected system 
dissipates in the projected phase ordering dynamics. To 
explain this dissipation thoroughly, we have to consider 
some degrees of freedom that were neglected in our anal- 
ysis, such as the motions of the "^-component along the 
vortex core, emissions of vortex rings, and density pulses 
to the extra dimension induced by vortex reconnection 
S4I. 
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FIG. 5: Decay of the projected vortex core density. Different 
marks show the plots of the density I20 for v = 0.80, 0.82, 
0.84, 0.86, 0.88, 0.90, and 0.92. Each plot obeys a power 
law t~ 1//2 when the structure factor is close to the scaling 
function, implying that the projected field behaves as a non- 
conserved field. The values of hu at time t = t v , after which 
I20 obeys the power law, are compared to the estimations 
of the defect nucleation rate in Fig. [4] For large v, the value 
l2u{t v ) increases unphysically and the plots fluctuate strongly 
because many ghost vortices remain because of a large amount 
of the ^2-component. 



